***************************** 
 * Paper: Terremotos y Crimen
 * Date: 9/Feb/2022
 *****************************

 ***************************
 * I. Preparacion de datos *
 ***************************

 * S1 Prepare Data.
 * Usaremos  8 meses antes y 15 despues del terremoto: 58,968 observaciones

 use "C:\Users\jose.balmori\Universidad Anahuac Dropbox\Jose Balmori\Paper 50 Terremotos\DataN\Base.dta"
 keep if ano==2017 | ano==2018  | ano==2019
 tab mes ano
 save "C:\Users\jose.balmori\Universidad Anahuac Dropbox\Jose Balmori\Paper 50 Terremotos\DataN\BaseProyecto.dta"

 * S2 Identificar municipios afectados por sismo de 2017

 use "C:\Users\jose.balmori\Universidad Anahuac Dropbox\Jose Balmori\Paper 50 Terremotos\DataN\BaseProyecto.dta"
 merge m:1 Clave using "C:\Users\jose.balmori\Universidad Anahuac Dropbox\Jose Balmori\Paper 50 Terremotos\DataN\Muncipios Terremoto 2017.dta"

   * Result                      Number of obs
   * -----------------------------------------
   * Not matched                        63,648
   *     from master                    63,648  (_merge==1)
   *     from using                          0  (_merge==2)

   * Matched                            24,804  (_merge==3)
   * -----------------------------------------


  keep if _merge==3 | _merge==1
  drop _merge
  
  *S3 Limpiar variables para identificacion. 
  
  recode earthquake (.=0)
  tab earthquake
  
  drop ent
  
  tostring Clave, gen(strClave)  /* Pasa la variable numerica a string */
  gen length = length(strClave)  /*  Auxiliar para manejar el tamaño de la Clave. Tamaño 4 y 5 */
  gen str8 var1=substr(strClave,1,1) if length==4
  gen str8 var2=substr(strClave,1,2) if length==5
  gen ents=var1 + var2
  drop var1 var2
  rename ents ent
  tab ent 
  
  destring, replace
  
  drop length strClave
  
  keep if ano==2017 | ano==2018 | ano==2019
  
  tab ano
  
  sort Clave
  quietly by Clave: gen dup= cond(_N==1,0,_n)
  tab dup /* Debes tener 24 observaciones por municipio */
  drop dup
  
  save "C:\Users\jose.balmori\Universidad Anahuac Dropbox\Jose Balmori\Paper 50 Terremotos\DataN\BaseProyecto.dta", replace
  
  *S4 Poner variables de interes en logaritmos.
  
  foreach i of varlist Homicidiost HomDolt HomCult Lesionest LesDolt LesCult Femt Abortot Secuestrot AbusoSexualt AcosoSexualt HostigamientoSexualt Violaciont ViolSimplet ViolEquit Robot RoboCasat RoboVeht RoboAutopartest RoboTrasnportistat RoboTranseuntet RoboTranseunteVPt RoboTranseunteEAt RoboTransportePubIndt RoboTransportePubColt RoboTransporteIndividualt RoboBancot RoboNegociot RoboGanadot RoboMaquinariat RoboOtrot Fraudet AbusoConft Extorsiont DanoPropt Despojot IPVt AsistFamt CorrMenort Narcot Amenazast Allanamientot Falsedadt Falsificaciont MedioAmbientet ServPubt Electoralest {
 	gen l`i'=log(1+`i')
 }
  
 * S5 Generar variable event. 
  gen even=.
  replace even=-8 if mes==1 & ano==2017 & earthquake==1
  replace even=-7 if mes==2 & ano==2017  & earthquake==1
  replace even=-6 if mes==3 & ano==2017   & earthquake==1
  replace even=-5 if mes==4 & ano==2017  & earthquake==1
  replace even=-4 if mes==5 & ano==2017  & earthquake==1
  replace even=-3 if mes==6 & ano==2017  & earthquake==1
  replace even=-2 if mes==7 & ano==2017  & earthquake==1
  replace even=-1 if mes==8 & ano==2017  & earthquake==1
  replace even=0 if mes==9 & ano==2017  & earthquake==1
  replace even=1 if mes==10 & ano==2017  & earthquake==1
  replace even=2 if mes==11 & ano==2017  & earthquake==1
  replace even=3 if mes==12 & ano==2017  & earthquake==1
  
  replace even=4 if mes==1 & ano==2018 & earthquake==1
  replace even=5 if mes==2 & ano==2018  & earthquake==1
  replace even=6 if mes==3 & ano==2018  & earthquake==1
  replace even=7 if mes==4 & ano==2018  & earthquake==1
  replace even=8 if mes==5 & ano==2018  & earthquake==1
  replace even=9 if mes==6 & ano==2018  & earthquake==1
  replace even=10 if mes==7 & ano==2018  & earthquake==1
  replace even=11 if mes==8 & ano==2018  & earthquake==1
  replace even=12 if mes==9 & ano==2018  & earthquake==1
  replace even=13 if mes==10 & ano==2018  & earthquake==1
  replace even=14 if mes==11 & ano==2018  & earthquake==1
  replace even=15 if mes==12 & ano==2018  & earthquake==1
  
  
  replace even=16 if mes==1 & ano==2019 & earthquake==1
  replace even=17 if mes==2 & ano==2019  & earthquake==1
  replace even=18 if mes==3 & ano==2019  & earthquake==1
  replace even=19 if mes==4 & ano==2019  & earthquake==1
  replace even=20 if mes==5 & ano==2019  & earthquake==1
  replace even=21 if mes==6 & ano==2019  & earthquake==1
  replace even=22 if mes==7 & ano==2019  & earthquake==1
  replace even=23 if mes==8 & ano==2019  & earthquake==1
  replace even=24 if mes==9 & ano==2019  & earthquake==1
  replace even=25 if mes==10 & ano==2019  & earthquake==1
  replace even=26 if mes==11 & ano==2019  & earthquake==1
  replace even=27 if mes==12 & ano==2019  & earthquake==1
  recode even (.=-1) /* Controles */
  tab even, m
  
   save "C:\Users\jose.balmori\Universidad Anahuac Dropbox\Jose Balmori\Paper 50 Terremotos\DataN\BaseProyecto.dta", replace
   
  * S6 Labels 
  
  sum lRoboCasat lRoboVeht lDanoPropt lAmenazast lLesionest lHomicidiost
  
   label var lAmenazast "Assault"
   label var lLesionest "Aggravated Assault"
   label var lHomicidiost "Homicides"

  label var lRoboCasat "Home Burglary"
  label var lRoboVeht  "Vehicle Theft"
  label var lDanoPropt "Vandalism"

 
 
 
  
 
  save "C:\Users\jose.balmori\Universidad Anahuac Dropbox\Jose Balmori\Paper 50 Terremotos\DataN\BaseProyecto.dta", replace
  
  * S7  Generate variable  After
  
  tab earthquake
  tab mes ano
  
  
  gen After=0
  replace After=1 if mes==9 & ano==2017  /*This is the month when the eartquke start */
  replace After=1 if mes==10 & ano==2017 
  replace After=1 if mes==11 & ano==2017  
  replace After=1 if mes==12 & ano==2017 
  
  replace After=1 if mes==1 & ano==2018 
  replace After=1 if mes==2 & ano==2018 
  replace After=1 if mes==3 & ano==2018 
  replace After=1 if mes==4 & ano==2018  
  replace After=1 if mes==5 & ano==2018  
  replace After=1 if mes==6 & ano==2018 
  replace After=1 if mes==7 & ano==2018  
  replace After=1 if mes==8 & ano==2018  
  replace After=1 if mes==9 & ano==2018    
  replace After=1 if mes==10 & ano==2018  
  replace After=1 if mes==11 & ano==2018  
  replace After=1 if mes==12 & ano==2018  
  
  
  replace After=1 if mes==1 & ano==2019 
  replace After=1 if mes==2 & ano==2019 %
  replace After=1 if mes==3 & ano==2019 
  replace After=1 if mes==4 & ano==2019  
  replace After=1 if mes==5 & ano==2019  
  replace After=1 if mes==6 & ano==2019 
  replace After=1 if mes==7 & ano==2019  
  replace After=1 if mes==8 & ano==2019  
  replace After=1 if mes==9 & ano==2019    
  replace After=1 if mes==10 & ano==2019  
  replace After=1 if mes==11 & ano==2019  
  replace After=1 if mes==12 & ano==2019  
  

  
  tab earthquake After
  
  save "C:\Users\jose.balmori\Universidad Anahuac Dropbox\Jose Balmori\Paper 50 Terremotos\DataN\BaseProyecto.dta", replace
  
  
   * S8 Año meses
   
  gen anomes = ym(ano, mes)
  format anomes %tm       
  save "C:\Users\jose.balmori\Universidad Anahuac Dropbox\Jose Balmori\Paper 50 Terremotos\DataN\BaseProyecto.dta", replace
  
  
  * S9 terremoto
  
  gen Terremoto= After*earthquake
  tab Terremoto
  label var Terremoto "Earthquake"
  
  
  
  ****************************
 * II. Tablas y Figuras     *
 ****************************
 
 ************************************
 * Table 1 Descriptive Statistics. 
 ************************************
 
 clear all
 use "C:\Users\jose.balmori\Universidad Anahuac Dropbox\Jose Balmori\Paper 50 Terremotos\DataN\BaseProyecto.dta"
 keep if ano==2017 | ano==2018
 
 *ssc install estout
 
 * Labels
 
  label var Amenazast "Assault"
   label var Lesionest "Aggravated Assault"
   label var Homicidiost "Homicides"

  label var RoboCasat "Residential Burglary"
  label var RoboVeht  "Vehicle Theft"
  label var DanoPropt "Vandalism"
 
 *Tables
 
 generate AfterA= After /* This is implemented in order to have the correct signs */
 recode AfterA (0=1) (1=0)
 
 sum  RoboVeht RoboCasat DanoPropt Amenazast Lesionest Homicidiost
 
 eststo contrA: estpost summarize RoboVeht RoboCasat DanoPropt Amenazast Lesionest Homicidiost  [aw=pob] if earthquake==0 & After==0
 
 eststo contrB: estpost summarize RoboVeht RoboCasat DanoPropt Amenazast Lesionest Homicidiost [aw=pob] if earthquake==0 & After==1
 
eststo controldif: estpost ttest RoboVeht RoboCasat DanoPropt Amenazast Lesionest Homicidiost if earthquake==0 , by(AfterA)   

 eststo treatA: estpost summarize RoboVeht RoboCasat DanoPropt Amenazast Lesionest Homicidiost [aw=pob] if earthquake==1 & After==0
 
 eststo treatB: estpost summarize RoboVeht RoboCasat DanoPropt Amenazast Lesionest Homicidiost [aw=pob] if earthquake==1 & After==1
 
 eststo treatdif: estpost ttest RoboVeht RoboCasat DanoPropt Amenazast Lesionest Homicidiost if earthquake==1, by(AfterA) 
 
 esttab contrA  contrB controldif treatA treatB treatdif  using Table1.tex, cell("mean(pattern(1 1 0 1 1 0) fmt(2)) b(star pattern(0 0 1 0 0 1)  fmt(2))") mgroups("Not Impacted" "Impacted",  pattern(0 0 0 1 0 0 ))  mtitle("Before" "After" "Difference" "Before" "After" "Difference") 
 
 *NOTA: Pero aca salen las diferencias con pesos que es lo que queremos
 
 foreach i of  varlist  RoboVeht RoboCasat DanoPropt Amenazast Lesionest Homicidiost {
 	reg `i'  After [aw=pob] if earthquake==0
 }
 
 foreach i of  varlist  RoboVeht RoboCasat DanoPropt Amenazast Lesionest Homicidiost {
 	reg `i'  After [aw=pob] if earthquake==1
 }
 
 
 
 *Tables: Sin pesos
 
 generate AfterA= After /* This is implemented in order to have the correct signs */
 recode AfterA (0=1) (1=0)
 
 sum  Amenazast Lesionest RoboVeht RoboCasat DanoPropt 
 
 eststo contrA: estpost summarize Amenazast Lesionest RoboVeht RoboCasat DanoPropt    if earthquake==0 & After==0
 
 eststo contrB: estpost summarize Amenazast Lesionest RoboVeht RoboCasat DanoPropt   if earthquake==0 & After==1
 
eststo controldif: estpost ttest Amenazast Lesionest RoboVeht RoboCasat DanoPropt  if earthquake==0 , by(AfterA)   

 eststo treatA: estpost summarize Amenazast Lesionest RoboVeht RoboCasat DanoPropt   if earthquake==1 & After==0
 
 eststo treatB: estpost summarize Amenazast Lesionest RoboVeht RoboCasat DanoPropt   if earthquake==1 & After==1
 
 eststo treatdif: estpost ttest Amenazast Lesionest RoboVeht RoboCasat DanoPropt  if earthquake==1, by(AfterA) 
 
 esttab contrA  contrB controldif treatA treatB treatdif  using Table1.tex, cell("mean(pattern(1 1 0 1 1 0) fmt(2)) b(star pattern(0 0 1 0 0 1)  fmt(2))") mgroups("Not Impacted" "Impacted",  pattern(0 0 0 1 0 0 ))  mtitle("Before" "After" "Difference" "Before" "After" "Difference") 
 
 *NOTA: Pero aca salen las diferencias con pesos que es lo que queremos
 
 foreach i of  varlist  Amenazast Lesionest RoboVeht RoboCasat DanoPropt {
 	reg `i'  After  if earthquake==0
 }
 
 foreach i of varlist  Amenazast Lesionest RoboVeht RoboCasat DanoPropt {
 	reg `i'  After  if earthquake==1
 }
 
 
 clear all
 
 ************************************
 * Table 2 Dif-in-Dif
 ************************************
 
 clear all
 use "C:\Users\jose.balmori\Universidad Anahuac Dropbox\Jose Balmori\Paper 50 Terremotos\DataN\BaseProyecto.dta"
 keep if ano==2017 | ano==2018
 

 
 sum lAmenazast lLesionest lHomicidiost lRoboVeht lRoboCasat lDanoPropt 
 
 foreach var of varlist lAmenazast lLesionest lHomicidiost lRoboVeht lRoboCasat lDanoPropt   {
 	qui reghdfe `var' Terremoto i.mes i.ano  [aw=pob], absorb(i.Clave) cluster(Clave)
est store D`var'
}
 
 esttab DlAmenazast DlLesionest DlHomicidiost DlRoboVeht DlRoboCasat DlDanoPropt ///
 using  "Tabla2A.tex", ///
 b(3) se(3) se label star(* 0.1 ** 0.05 *** 0.01) ///
 title(Difference-in-Differences Specification:   \label{lpmpsycho}) ///
 nomtitles r2(2) obslast 
 
 
 
 
 
 
 
 ************************************
 * Fig III  Event Study
 ************************************
 
  clear all
  use "C:\Users\jose.balmori\Universidad Anahuac Dropbox\Jose Balmori\Paper 50 Terremotos\DataN\BaseProyecto.dta"
  keep if ano==2017 | ano==2018
  
  cd "C:\Users\jose.balmori\Universidad Anahuac Dropbox\Jose Balmori\Paper 50 Terremotos\OutputN"
  label var lRoboCasat "Residential Burglary"
  
  grstyle init
  grstyle set lpattern short-dash  solid     /*patron de las dos lineas */
  grstyle set linewidth    medthick  thick /*Grosor de las dos lineas */
  grstyle set plain  , horizontal compact  minor dotted
  grstyle set color  "56 108 176"    "56 108 176"    
  grstyle set symbol 
 
  local dep   lAmenazast lLesionest lRoboVeht lRoboCasat lDanoPropt
 
 foreach v of local dep {
 local label : variable label `v'  
 eventdd `v' i.mes i.ano  [aw=pob],  hdfe absorb(Clave) timevar(even) cluster(i.Clave) ci(rline) ///
			graph_op( xlabel(-8(1)15)  /*change month */ xtitle("Months since the earthquake") ytitle("Change in Crimes per 100,000 Inhabitants")  name(IPVt , replace) title("`label'", pos(12)) recast(line) )  
 graph save `v', replace
 }
 
 * Cambio Manual
 * Ir a Graph Editor -> plotregion1 (1 click)
 * Paso 1: Ir a References lines -> line y -> Color: Negro
 * paso 2: Ir a plot2 ->Connected -> Tamaño: Mediano -> Simbolo (Circulo)
 

grc1leg  lAmenazast.gph lLesionest.gph lRoboVeht.gph lRoboCasat.gph lDanoPropt.gph       ///
, iscale(.5) col(2)    imargin(0 0 0 0 ) name(g7, replace)  
graph display g7, xsize(6) ysize(8)
graph save F1, replace
graph export "C:\Users\jose.balmori\Universidad Anahuac Dropbox\Jose Balmori\Paper 50 Terremotos\OutputN\Fig1.png", as(png) width(1000) replace


 ****************************
 * III. Robustness     *
 ****************************
 
 
 * Tabla A1.A Functional Form (lineal)
 
  clear all
  use "C:\Users\jose.balmori\Universidad Anahuac Dropbox\Jose Balmori\Paper 50 Terremotos\DataN\BaseProyecto.dta"
  keep if ano==2017 | ano==2018
    
 
 foreach var of varlist Amenazast Lesionest Homicidiost RoboVeht RoboCasat DanoPropt    {
 	qui reghdfe `var' Terremoto i.mes i.ano [aw=pob], absorb(i.Clave) cluster(Clave)
est store D`var'
}
 
 esttab DAmenazast DLesionest DHomicidiost DRoboVeht DRoboCasat DDanoPropt  ///
 using  "TablaR1a.tex", ///
 b(3) se(3) se label star(* 0.1 ** 0.05 *** 0.01) ///
 title(Difference-in-Differences Specification:   \label{lpmpsycho}) ///
 nomtitles r2(2) obslast 
 
 
  
 * Tabla A1.B Without Weights
 
  clear all
  use "C:\Users\jose.balmori\Universidad Anahuac Dropbox\Jose Balmori\Paper 50 Terremotos\DataN\BaseProyecto.dta"
  keep if ano==2017 | ano==2018
  
 
 foreach var of varlist lAmenazast lLesionest lHomicidiost lRoboVeht lRoboCasat lDanoPropt    {
 	qui reghdfe `var' Terremoto i.mes i.ano, absorb(i.Clave) cluster(Clave)
est store D`var'
}
 
 esttab DlAmenazast DlLesionest DlHomicidiost DlRoboVeht DlRoboCasat DlDanoPropt  ///
 using  "TablaR1b.tex", ///
 b(3) se(3) se label star(* 0.1 ** 0.05 *** 0.01) ///
 title(Difference-in-Differences Specification:   \label{lpmpsycho}) ///
 nomtitles r2(2) obslast 
 
  
 
 
 * Tabla A1.C Without Mexico City
 
  clear all
  use "C:\Users\jose.balmori\Universidad Anahuac Dropbox\Jose Balmori\Paper 50 Terremotos\DataN\BaseProyecto.dta"
  keep if ano==2017 | ano==2018
  
 
 foreach var of varlist lAmenazast lLesionest lHomicidiost lRoboVeht lRoboCasat lDanoPropt    {
 	qui reghdfe `var' Terremoto i.mes i.ano [aw=pob] if ent~=9, absorb(i.Clave) cluster(Clave)
est store D`var'
}
 
 esttab DlAmenazast DlLesionest DlHomicidiost DlRoboVeht DlRoboCasat DlDanoPropt   ///
 using  "TablaR1c.tex", ///
 b(3) se(3) se label star(* 0.1 ** 0.05 *** 0.01) ///
 title(Difference-in-Differences Specification:   \label{lpmpsycho}) ///
 nomtitles r2(2) obslast 
 
 * Tabla A1.D Multiple Hypothesis testing
 

 
  clear all
  use "C:\Users\jose.balmori\Universidad Anahuac Dropbox\Jose Balmori\Paper 50 Terremotos\DataN\BaseProyecto.dta"
  keep if ano==2017 | ano==2018
  areg lAmenazast Terremoto i.mes i.ano [aw=pob], absorb(Clave) cluster(Clave)
  regsave, pval
  keep var pval
  keep if pval!=.
  save "C:\Users\jose.balmori\Universidad Anahuac Dropbox\Jose Balmori\Paper 50 Terremotos\OutputN\pval1.dta"
  
  
   clear all
  use "C:\Users\jose.balmori\Universidad Anahuac Dropbox\Jose Balmori\Paper 50 Terremotos\DataN\BaseProyecto.dta"
  keep if ano==2017 | ano==2018
  areg lLesionest Terremoto i.mes i.ano [aw=pob], absorb(Clave) cluster(Clave)
  regsave, pval
  keep var pval
  keep if pval!=.
  save "C:\Users\jose.balmori\Universidad Anahuac Dropbox\Jose Balmori\Paper 50 Terremotos\OutputN\pval2.dta"
  
  clear all
  use "C:\Users\jose.balmori\Universidad Anahuac Dropbox\Jose Balmori\Paper 50 Terremotos\DataN\BaseProyecto.dta"
  keep if ano==2017 | ano==2018 
  areg lHomicidiost Terremoto i.mes i.ano [aw=pob], absorb(Clave) cluster(Clave)
  regsave, pval
  keep var pval
  keep if pval!=.
  save "C:\Users\jose.balmori\Universidad Anahuac Dropbox\Jose Balmori\Paper 50 Terremotos\OutputN\pval3.dta"
  
  
  
   clear all
  use "C:\Users\jose.balmori\Universidad Anahuac Dropbox\Jose Balmori\Paper 50 Terremotos\DataN\BaseProyecto.dta"
  keep if ano==2017 | ano==2018
  areg lRoboVeht Terremoto i.mes i.ano [aw=pob], absorb(Clave) cluster(Clave)
  regsave, pval
  keep var pval
  keep if pval!=.
  save "C:\Users\jose.balmori\Universidad Anahuac Dropbox\Jose Balmori\Paper 50 Terremotos\OutputN\pval4.dta"
  
  
   clear all
  use "C:\Users\jose.balmori\Universidad Anahuac Dropbox\Jose Balmori\Paper 50 Terremotos\DataN\BaseProyecto.dta"
  keep if ano==2017 | ano==2018
  areg lRoboCasat  Terremoto i.mes i.ano [aw=pob], absorb(Clave) cluster(Clave)
  regsave, pval
  keep var pval
  keep if pval!=.
  save "C:\Users\jose.balmori\Universidad Anahuac Dropbox\Jose Balmori\Paper 50 Terremotos\OutputN\pval5.dta"
  
   clear all
  use "C:\Users\jose.balmori\Universidad Anahuac Dropbox\Jose Balmori\Paper 50 Terremotos\DataN\BaseProyecto.dta"
  keep if ano==2017 | ano==2018
  areg lDanoPropt  Terremoto i.mes i.ano [aw=pob], absorb(Clave) cluster(Clave)
  regsave, pval
  keep var pval
  keep if pval!=.
  save "C:\Users\jose.balmori\Universidad Anahuac Dropbox\Jose Balmori\Paper 50 Terremotos\OutputN\pval6.dta"
  
 
 
  use "C:\Users\jose.balmori\Universidad Anahuac Dropbox\Jose Balmori\Paper 50 Terremotos\OutputN\pval1.dta"
  append using "C:\Users\jose.balmori\Universidad Anahuac Dropbox\Jose Balmori\Paper 50 Terremotos\OutputN\pval2.dta"
  append using "C:\Users\jose.balmori\Universidad Anahuac Dropbox\Jose Balmori\Paper 50 Terremotos\OutputN\pval3.dta"
  append using "C:\Users\jose.balmori\Universidad Anahuac Dropbox\Jose Balmori\Paper 50 Terremotos\OutputN\pval4.dta"
  append using "C:\Users\jose.balmori\Universidad Anahuac Dropbox\Jose Balmori\Paper 50 Terremotos\OutputN\pval5.dta"
  append using "C:\Users\jose.balmori\Universidad Anahuac Dropbox\Jose Balmori\Paper 50 Terremotos\OutputN\pval6.dta"
  save "C:\Users\jose.balmori\Universidad Anahuac Dropbox\Jose Balmori\Paper 50 Terremotos\OutputN\pval.dta"
 
  cd "C:\Users\jose.balmori\Universidad Anahuac Dropbox\Jose Balmori\Paper 50 Terremotos\OutputN\"
  
  * Open: 5. fdr_sharpended_qvalues.do and run since line 37
  * The column bky06_qval provides the q values. 
  
  save "C:\Users\jose.balmori\Universidad Anahuac Dropbox\Jose Balmori\Paper 50 Terremotos\OutputN\pval.dta", replace
  
  
  * Tabla A1.E Different control group
 
  clear all
  use "C:\Users\jose.balmori\Universidad Anahuac Dropbox\Jose Balmori\Paper 50 Terremotos\DataN\BaseProyecto.dta"
  keep if ano==2017 | ano==2018
  
  gen experimento2=0 
  replace experimento2=1 if ent==20 | ent==7 | ent==30 |ent==9 | ent==17 | ent==21 | ent==15 | ent==29 | ent==12   ///
  | ent==14  | ent==6  | ent==18  | ent==1  | ent==11  | ent==22  | ent==13  | ent==27 | ent==31  | ent==4  | ent==23
   tab experimento2
  
 
 foreach var of varlist lAmenazast lLesionest lHomicidiost lRoboVeht lRoboCasat lDanoPropt    {
 	qui reghdfe `var' Terremoto i.mes i.ano [aw=pob] if experimento2==1, absorb(i.Clave) cluster(Clave)
est store D`var'
}
 
 esttab DlAmenazast DlLesionest DlHomicidiost DlRoboVeht DlRoboCasat DlDanoPropt  ///
 using  "TablaR1e.tex", ///
 b(3) se(3) se label star(* 0.1 ** 0.05 *** 0.01) ///
 title(Difference-in-Differences Specification:   \label{lpmpsycho}) ///
 nomtitles r2(2) obslast 
 
  * Tabla A1.F  Incluyendo como control municipios en estados no afectados
  
  
  clear all
  use "C:\Users\jose.balmori\Universidad Anahuac Dropbox\Jose Balmori\Paper 50 Terremotos\DataN\BaseProyecto.dta"
  keep if ano==2017 | ano==2018
 
 * generar variable de municipiuos a excluir
 
  tab ent if earthquake==1
  
  gen enteartquake=0
  replace enteartquake=1 if ent==7|ent==9 |ent==12 |ent==15 |ent==17 |ent==20 |ent==21 |ent==29 |ent==30 
  tab enteartquake
  
  gen  incluir=1
 replace incluir=0 if  enteartquake==1 & earthquake==0 /*Estos son los municipios de estados adectados pero sin terremoto*/
 tab incluir
  
  foreach var of varlist lAmenazast lLesionest lHomicidiost lRoboVeht lRoboCasat lDanoPropt    {
 	qui reghdfe `var' Terremoto i.mes i.ano [aw=pob] if incluir==1, absorb(i.Clave) cluster(Clave)
est store D`var'
}
 
 esttab DlAmenazast DlLesionest DlHomicidiost DlRoboVeht DlRoboCasat DlDanoPropt  ///
 using  "TablaR1f.tex", ///
 b(3) se(3) se label star(* 0.1 ** 0.05 *** 0.01) ///
 title(Difference-in-Differences Specification:   \label{lpmpsycho}) ///
 nomtitles r2(2) obslast 
  
  
  * Tabla A1.G  Espandiendo el tiempo
  
  
 clear all
 use "C:\Users\jose.balmori\Universidad Anahuac Dropbox\Jose Balmori\Paper 50 Terremotos\DataN\BaseProyecto.dta"
 tab ano
 tab Terremoto
 

 
 sum lAmenazast lLesionest lHomicidiost lRoboVeht lRoboCasat lDanoPropt 
 
 foreach var of varlist lAmenazast lLesionest lHomicidiost lRoboVeht lRoboCasat lDanoPropt   {
 	qui reghdfe `var' Terremoto i.mes i.ano  [aw=pob], absorb(i.Clave) cluster(Clave)
est store D`var'
}
 
 esttab DlAmenazast DlLesionest DlHomicidiost DlRoboVeht DlRoboCasat DlDanoPropt ///
 using  "TablaR1g.tex", ///
 b(3) se(3) se label star(* 0.1 ** 0.05 *** 0.01) ///
 title(Difference-in-Differences Specification:   \label{lpmpsycho}) ///
 nomtitles r2(2) obslast 
 
 